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Abstract 

In high-energy processes near the endpoint, there emerge new contributions associated with 
spectator interactions. Away from the endpoint region, these new contributions are suppressed 
compared to the leading contribution, but the leading contribution becomes suppressed as we ap- 
proach the endpoint and the new contributions become comparable. We present how the new 
contributions scale as we reach the endpoint and show that they are comparable to the suppressed 
leading contributions in deep inelastic scattering by employing a power counting analysis. The 

*y-N hadronic tensor in deep inelastic scattering is shown to factorize including the spectator interac- 

tions, and it can be expressed in terms of the lightcone distribution amplitudes of initial hadrons. 

i We also consider the contribution of the spectator contributions in Drell-Yan processes. Here the 

spectator interactions are suppressed compared to double parton annihilation according to the 
power counting. 
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I. INTRODUCTION 

The standard results of high-energy scattering processes based on the operator product 
expansion are consistent and work very well, but they are expected to be modified near 
the endpoint where the Bjorken variable x approaches 1. Since the available phase space is 
restricted near the endpoint region, peculiar physical results arise and there has been a lot 
of theoretical interest in the endpoint behavior of high-energy scattering processes. 

The kinematic peculiarity near the endpoint x ~ 1 — A/Q manifests two features which 
do not show up away from the endpoint region, where Q is a large scale and A is the 
typical QCD scale for hadron masses. Firstly, the soft Wilson lines accompanied by collinear 
particles do not cancel completely, and the remnant is combined to produce soft functions. 
Extracting the soft part is crucial in factorization proof. In deep inelastic scattering (DIS) 
near the endpoint, since the invariant mass of the final-state particles is p 2 x ~ <5 2 (1 — x), 
spectator particles after hard scattering can be either soft or collinear to the final-state jets 
leaving no particles in the beam direction. In Drell-Yan (DY) process near the endpoint 
r = Q 2 /s — ¥ 1, there can be only soft final-state particles except a lepton pair due to the 
kinematic constraint. Comparing these two processes near the endpoint, the configurations 
of soft particles in the final states are different, causing different types of soft interactions, 
and the factorization proof near the endpoint is affected significantly by the soft parts. 

Secondly, the contribution of spectator partons to the scattering cross section, which is 
subleading away from the endpoint region, is not negligible near the endpoint and it should 
be included in the scattering cross section. It is not because the spectator contributions are 
enhanced, but because the leading contribution is suppressed near the endpoint to become 
of the same order as the spectator contribution. The proof that the spectator interaction 
becomes also important and the factorization property including the spectator interactions 
are the main theme of this paper. 

The momentum of an energetic hadron in the lightlike n-direction can be decomposed 
into 

l? = n-p—+rt_ + n-p— = 0(Q) + 0(A) + 0(A 2 /Q), (1) 

where the lightcone vectors n M and n M satisfy n 2 = n 2 = and n ■ n — 2. The hadron is 
constrained to be on the mass shell p 2 ~ A 2 , so are the partons constituting the hadron, such 
that a scattering process can be described in terms of the parton distribution functions (PDF) 
as the probability distribution. However, these constraints give rise to special kinematic 
situation near the endpoint. Since the active parton undergoing hard scattering carries 
most of the energy inside the hadron, the n-component of the momentum for the spectator 
partons is of order A. These spectator partons can have momenta satisfying the relative 
scaling to be n-collinear, but they cannot be on the mass shell. If the spectator partons 
become soft with all the momentum components of order A, they can be on the mass shell. 
But the total momentum of the hadron, being the sum of a collinear and a soft momenta, 
becomes of order P 2 ~ QA, which is far off mass shell. Therefore in order to be consistent 
with the constraints of the on-shellness at the partonic and at the hadronic levels, and 
the kinematic constraint in the endpoint region, the initial spectator quarks are energetic, 



n-collinear, and undergo a large momentum transfer inside the hadron of order Q 2 or QA. 

As a result, near the endpoint region, the spectator particles which are initially n-collinear 
become either n-collinear or soft after the large momentum transfer of order Q 2 or QA. This 
momentum transfer is related not to the hard scattering, but to the spectator interaction in 
the initial hadron. This necessitates the spectator interaction with a large momentum trans- 
fer in the scattering process near the endpoint in order to reflect the kinematic restrictions 
consistently. 

As we will see later, it is the main reason for the suppression of the conventional scatter- 
ing cross section [j near the endpoint, which becomes comparable to the contribution of the 
spectator interactions. In the standard region 1 — x ~ 0(1), the spectator contribution with 
a large momentum transfer is suppressed by A 2 /Q 2 compared to the leading conventional 
contribution, thus can be safely neglected. All-order factorization analyses (in a s ) were pre- 
sented in Refs. [U |2] near the endpoint region for Drell-Yan processes, and the subleading 
contributions suppressed by powers of A/Q from the final-state interactions via the sublead- 
ing final-state jet functions were analyzed in Refs. [?H9]. However, the issue of the spectator 
contribution has not been addressed in the limit x — > 1 in previous literature. Careful power 
counting indicates that the leading contribution obtained away from the endpoint region 
experiences severe suppression such that it is comparable to the spectator contributions as 
x goes to 1. 

DIS in the endpoint region has been so far conventionally described by the following 
schematic factorization formula pfl E] 

F,(Q 2 , x) ~ H(Q 2 , fM) ■ J(Q 2 (1 -x),im)® f i/H (x, /i), (2) 

where F\ is the conventional structure function in the endpoint, H is a hard function, and 
fi/H is a PDF. J represents the final-state jet function integrating out the degrees of freedom 
of order Q 2 (l — x). And '(g)' denotes the convolution of the jet function with the PDF. In 
the framework of soft-collinear effective theory (SCET) [3H5], this factorization formula 
has been revisited and confirmed without considering the spectator interactions [2HSJ- If 
the spectator contribution should be included near the endpoint as discussed above, the 
conventional leading contribution of Eq. (J2J) is to be modified including this contribution, 
too. The PDF includes both the collinear part in the beam and the soft part. The collinear 
part can be described by the lightcone distribution amplitudes (LCDA) for the initial hadron, 
and the soft part includes the final-state soft spectator quarks, which modifies the structure 
of the PDF. 

This mechanism also affects the longitudinal structure function Fl- The dominant spec- 
tator contribution to Fl comes from the subleading corrections to the current operator re- 
sponsible for spectator interactions, and remarkably it becomes comparable to Fi, since F\ is 
suppressed near the endpoint. The Callan-Gross relation states that Fl = — Fi + F 2 Q 2 / (4:X 2 ) 
vanishes at leading order in 1/Q, but it does not have to hold at subleading order we con- 
sider here. If we consider the subleading jet function related to the final-state particles alone 



1 We mean the 'conventional scattering cross section' by the scattering cross section neglecting spectator 
partons. 



without the spectator contribution, it is shown that the contribution to Fl is suppressed by 
A/Q [7H§] compared to F x . This arises from the subleading jet function by integrating out 
the degrees of freedom of order p 2 ~ QA in the final state. However the new contributions 
which will be considered here turn out to be dominant, compared to the contribution to Fl 
from subleading jet functions without the spectator interaction. 

In Drell-Yan processes, the spectator particles can be in the original direction of the initial 
hadron as in DIS away from the endpoint region. Near the endpoint, since the invariant 
mass of the final-state hadrons is of order A 2 , there can be only soft particles. This is in 
contrast to DIS, since the spectator particles are either n-collinear (collinear to the final-state 
energetic collinear particle) or soft in DIS near the endpoint. In DIS, the case with final 
n-collinear particles corresponds to the endpoint limit of the conventional approach, and can 
be compared to the new contribution with the spectator interaction. But there is no such 
analog of the case with n-collinear particles in DY processes. However, the situation gets 
more drastic since we should also consider the double parton annihilation in DY process, for 
it is less suppressed than the spectator interaction as far as the power counting is concerned. 

There is one hadronic scalar function to describe DY processes. The soft part differs 
from that in DIS, hence needs some modification or a different definition in the PDF. Away 
from the endpoint region, the soft part cancels, and the PDF consists of the matrix elements 
of collinear operators. It enables us to use universal PDFs independent of the scattering 
processes. That is, if we obtain or define the PDF in DIS, it can be used in DY processes. 
Near the endpoint, the soft part does not cancel, and it should be included in the definition 
of the PDF. If the PDF defined in DIS is to be employed in DY processes, there should be 
some modification which incorporates the difference of the soft parts in the two processes. 

In this paper, we consider the new contributions arising from spectator interactions in 
DIS and DY processes. The power counting is performed systematically and it is shown 
that the size of the new contributions is comparable to the standard contribution near the 
endpoint and the factorization property is considered. In Section [IT] we perform the power 
counting analysis in DIS in the large x limit and show how the spectators engage in the 



scattering process. In Sec. Ill, we show the factorization property for the new contributions 



in DIS. We employ two-step matching to prove the factorization explicitly. In Sec. IV we 
present the power counting analysis for Drell-Yan processes, including the double parton 
annihilation. In Sec. |Vl we give a conclusion. 



II. POWER COUNTING IN DIS AS x -> 1 

A systematic power counting can be applied to study the suppression of the scattering 
cross section near the endpoint. Let us illustrate how the power counting is performed 
in DIS first. The momentum of the final states is given by px = Q + P, where q is the 
momentum transfer from the leptonic system and P is the momentum of the initial hadron. 
The invariant mass of the final states is given by 

Px= {1 ^ 1 Q 2 + m] I , (3) 

x 



where x = —q 2 /2P ■ q = Q 2 /2P ■ q and inn is the initial hadron mass. We now choose the 
Breit frame in which q^ = Q(n M — n M )/2 and the initial hadron is described as an n-collinear 
particle. The invariant mass of the final state varies as x changes. Away from the endpoint, 
Px ~ Q 2 1 an d this represents general hard scattering processes. As x — > 1, the invariant 
mass gets smaller, and the limit is classified into two regions. The first is the resonance 
region where 1 — x ~ A 2 /Q 2 with p 2 x ~ A 2 , in which only n-collinear hadrons are allowed 
kinematically in the final state. And the second is the endpoint region 1 — x ~ A/Q with 
Px ~ Q^i m which there can be n-collinear jets and soft hadrons. 

In both regions there are no n-collinear final-state particles, while the initial-state par- 
tons are n-collinear particles. Therefore the spectator particles have to interact with large 
momentum transfer to become either soft or n-collinear, and we have to include all the 
interactions of the initial partons. In the resonance region all the spectators undergo hard 
interactions with the momentum transfer of order Q 2 , and then the spectator particles, which 
are initially n-collinear, are converted into n-collinear particles to make p\ ~ A 2 . In the 
endpoint region we have two possibilities: First, an n-collinear spectator inside the initial 
hadron can be n-collinear undergoing hard interactions as in the resonance region. But here 
the offshellness of the final state is allowed to be of order QA, much larger than the resonant 
case. Secondly, a spectator loses most of its energy to the active parton and becomes a soft 
particle. This energy transfer between the active parton and the spectator is hard-collinear 
in the n direction. Its offshellness is of order QA, which is the typical offshellness of the 
final-state jet in the endpoint region. 

Now we can perform the power counting of the hadronic tensor for inclusive DIS, which 
is defined as 

W» v = E / '^■ z (H\J^(z)\X)(X\r(0)\H), (4) 

where H is the initial hadron, J^ is an electromagnetic current, and the summation in- 
cludes the phase space of the final-state particles. The power counting on the volume 
d 4 z depends on how much phase space is available. For power counting on the remaining 
part J2x{H\J^\X)(X\J u \H), we divide it into three parts; the initial state, the amplitude 
squared, and the final state contributions. Here we focus on the amplitudes at tree level, 
but the result on the power counting by rj ~ A/Q can be easily extended to loop corrections 
because no loop contribution can enhance the amplitude by inverse powers of 77. 

In the standard region where 1 — x ~ 0{1), d 4 z covers the full phase space, hence power- 
counted as 1/Q 4 - The initial-state part is schematically written as |(0|\l/ n |if)| 2 . Here \l/ ra is 
an n-collinear quark and scales as Q 3 ^ 2 i] with 77 ~ A/Q, and the collinear state \H) scales 
as 1/A ~ 1/(Qtj). Therefore the initial-state part yields the factor Q. The final state 
contains J d 4 p5(p 2 )ft, where ft comes from the spinor sum of the final state. Because the 
final state carries the hard momentum in the standard region, the power counting states 
that J d 4 p5{p 2 )ft ~ Q 4 ■ (1/Q 2 ) • Q = Q 3 . Also the amplitude squared is simply 0(1). In the 
standard region, spectator contributions do no change the power counting since they are of 
order 1 or give higher powers of 77. Therefore the overall power counting for the structure 
function yields 0{1). 

As explained above, the spectator contribution should be included near the endpoint 



region. It is also important how many particles there are in the leading Fock space of 
the initial hadron H. In the case of a pion, there are qq in the leading Fock space and 
qqq for a proton. And the power counting on the structure functions for a pion and a 
proton is different. Because all the partons are involved in the scattering process, the time- 
ordered products of the electromagnetic current and the interaction Lagrangians including 
all the spectators should be taken into account in the hadronic tensor W^ v . For the power 
counting of the initial-state contributions, we consider | (Ol^^^I/^lTr) | 2 for an initial-state pion 
and |(0|\I/ n \I/ n \I/ ra |p)| 2 for a proton neglecting irrelevant Lorentz structure and color factors. 
From our power counting rule, these yield the factors Q 2 A 2 and Q 3 A 4 respectively, and the 
structure function for the proton is more suppressed than the structure function for the pion 
near the endpoint. 

In the resonance region the phase space is severely constrained and the invariant mass 
of the final states becomes p\ ~ Q 2 (l — x) + m 2 H ~ A 2 . The momentum px flows between 
the two points and z in the hadronic tensor, and it implies that the volume d A z is counted 
as 1/A 4 . Some examples of DIS near the endpoint for an initial pion and a proton are 
shown in Fig. [1} The momentum transfer between the active and the spectator quark is 
hard (p 2 ~ Q 2 ), hence the amplitudes for an initial pion and a proton scale as l/Q 3 and 
l/Q 6 respectively. Each quark field in the final state is power counted as / d i p c 5(p 2 )p' c ~ 
A 4 • (1/A 2 ) ■ Q = QA 2 , where p c represents collinear momentum with the offshellness of order 
A 2 . Combining all the factors, the power counting of the hadronic tensor is given as 

W^ ~ M 2 -I-F-V, 

£)».q.a'. W a')»-^~^~i-. *** = ,, (g) 

(^) 2 • Q 3 A 4 • (QA 2 ) 3 • ji ~ f ~ (1 - 4 3 for H=p, 

where Ai 2 denotes the amplitude squared, 1(F) is the initial (final) state, and V indicates 



<f = Q(W - n")/2 




n, soft 




FIG. 1. Specific examples of DIS processes (a) for a pion and (b) for a proton in the initial state 
near the endpoint. In the resonance region 1 — x ~ 0(rj 2 ), all the final-state quarks are n-collinear 
(p 2 ~ A 2 ). In the endpoint region 1 — x ~ 0(r)), the spectator quarks in the final state can be 
either n-hard-collinear (j? hc ~ QA) or soft (p 2 ~ A 2 ). 
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TABLE I. The scaling behavior of the hadronic tensor W^ v in the endpoint region (1 — x ~ A/Q). 
H is an initial hadron, ^>n,hc represents an n-hard-collinear quark, and q s is a soft quark. 

the volume d 4 z. This is consistent with the previous power counting in the resonance region 

In the endpoint region p 2 x scales as Q 2 (l — x) ~ QA and d A z is counted as 1/(QA) 2 . 
The spectator quarks in the final state can be either hard-collinear (p 2 c ~ QA) or soft 
(p 2 ~ A 2 ), while the active parton in the final state is kept to be hard-collinear for the 
maximal scaling. We estimate the power counting of amplitudes from Fig. [TJ For an initial 
pion the amplitude is power counted as either 1/Q 3 (hard momentum transfer) or 1/(Q 2 A) 
(hard-collinear momentum transfer). For an initial proton the amplitude is estimated to 
be of order 1/Q 6 (two hard-collinear spectators), 1/(Q 5 A) (one hard-collinear and one soft 
spectators), and 1/(Q 4 A 2 ) (two soft spectators). The hard-collinear final state is maximally 
power counted as / d A ph c 5{p 2 hc )f hc ~ Q 2 A 2 ■ (1/Q A) ■ Q = Q 2 A and the final soft state scales 
as / d 4l p s S(p 2 )p s ~ A 4 • (1/A 2 ) • A = A 3 . The final results of the scaling behavior of W^ v 
are summarized in Table, [ij Near the endpoint, the hadronic tensor scales as (1 - x) 2 for 
an initial pion, and (1 



x\ 



for an initial proton for all the possible final states. The point 
is that the suppression of the hadronic tensor near the endpoint is the same for the final 
n-collinear and soft particles, and depends only on the type of the initial hadrons. 



III. FACTORIZATION ANALYSIS OF DIS NEAR THE ENDPOINT 

In this section we analyze the factorization of DIS near the endpoint in SCET. For 
simplicity we consider DIS with an initial pion rather than with a proton. But the extension 
to the initial proton is straightforward. The general tensor structure of W^ u for DIS in the 
Breit frame can be written as 



W^ v 



sT + 



q H q u 



F 1 + (P 11 



P-q 



P" 



P-q , 
qz 



-gFF l + v*v v F L , (6) 



where g^/ = g^ u — (n^n u + n M n I/ )/2, and v M = (n M + n fl )/2. The longitudinal structure 
function Fl is defined as Fl = —F\ + F2Q 2 '/(Ax 2 ). Away from the endpoint region Fl is 
suppressed compared to F x . But as x goes to 1, F 1 is suppressed and F L becomes comparable 



to F\. Both structure functions are influenced by the spectator interaction in the endpoint 
region 1 — x ~ A/Q. 



A. Contribution from hard gluon exchange 

For the hard gluon exchange as shown in Fig. [TJ we obtain the local SCET operators 
with n and n-collinear quark fields by integrating out hard gluons. For an initial pion, these 
operators are obtained from Fig. [I] (a) along with the hard gluon exchange between the 
outgoing active quark and the spectator quark. After matching these contributions onto 
SCET, the electromagnetic current J M = q^q is expressed in terms of a convolution as 



v^C H ®O h = f m / dudw C H (u, w, Q)O h (u, w) 
Jo 

where Oh(u,w) is given by 



(7) 



H (u,w) 



Q 



-^5(w- n ^)YhlT a Yj( 



u 



n-V 



•qj . -qj y^ ± T a Y-'i!- 



(8) 



Here we take the active quark as a quark and the spectator quark as an antiquark. The 
SCET collinear field \P n (n) = W„(s)^n(n) is collinear-gauge invariant, where W n (n) is an n(n)- 
collinear Wilson line [I] . The variables u and w are the momentum fractions of the active 
quark before and after the hard scattering respectively, and n-V (n-lZ) is the label momentum 
of the n(n)-collinear field. Ch is the Wilson coefficient for the hard gluon exchange, and at 
tree level it is given by C H = 8Tca s /(uw) with u = 1 — u and w = I — w. 

We have redefined the collinear quark fields to decouple soft interactions as \I/ n (n) — > 
Y n (n)^ n (n) (annihilated quark) and \P n (n) — > Y n (n)^n(n) (created antiquark) in Eq. (p), where 
the soft Wilson lines are defined as (5J [H] 



Y n (x) = Pexp 
Y n (x) = Pexp 



ig I dsn ■ A s (ns) 

oo 



Y n (x) = Pexp —ig / dsn ■ A s (ns 



>fj 



dsn ■ A s (ns] 



(9) 



Yn{x) =Pexp — ig \ dsn ■ A s (ns) , (10) 



where P and P represent path-ordering and antipath-ordering respectively. Though there 
is only an octet four-quark operator at tree level, there can be singlet operators at higher 
order in a s , and we can take the appropriate color projection for a color-singlet pion. Note 
that the result in Eq. ([7]) is proportional to v^. Thus the hard-gluon exchange for the pion 
contributes to the longitudinal structure function Fl. 

If we take the matrix element of Eq. ([7]) between the initial pion and the final state X, 
the n-collinear part can be expressed in terms of the pion light-cone distribution amplitude 
(LCDA) [12] because there is no outgoing final n-collinear particle. Using the expression for 
the leading-twist LCDA in SCET 



<0| S(u 



n-V 
n-Pn 
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K(Pir)) 



1 f~ 5 



ab .^ 



(W^ 1 



U 



(11) 



we have 



u 



n 



n\i 



{X\C H ®0 H \it)=i^ 2 I dudwH(u,w,Q 2 )Mu)(Xn\^nS(w- ^7^)^75^10), (12) 
Q z Jo Q 2 

where (j) n is the leading-twist LCDA for the pion, N is the number of colors, and H is the 
hard factor given by CpCjj /2N at tree level. This expression can be generalized to include 



higher-order a s corrections. In Eq. (12) we put n ■ p w = Q neglecting (9(1 — x). The soft 



matrix element in Eq. ( 12 ) can be explicitly given as 



(X\C h ®O h \tx) = —^ £ dudwH(u,w,Q 2 )<PM 



Wilson lines in Eq. (|8|) cancel since the pion is a color singlet. By inserting Eq. (11), the 

(13) 

x (X^n8(u - ^ r )7^7 57±Q l^ t T a y n F n tT a F^ 7r |0), 

from which the cancellation of the soft Wilson lines can be clearly seen. 

From Eq. (El), the contribution of the hard-gluon exchange to Fl can be written as 

Ff (Q 2 ,x) = (2tt) 3 (^) 2 j 1 du'dw'dudwH*{u\w' ,Q 2 )^{u')H{u,w,Q 2 )(t)M (14) 



Q' 



o 



x J26(q + p« -p x )(0\*J(w' - ^)^ T5 * s |*n)<^|<M(™ - T ^)^75^|0), 



Q 



Q 



where the final-state jet function J^ is defined as 

Q 2 J T^4-^{ ( l + V^-Px)j"{w,w\p 2 x ) 

= E *(*+** - P x)(0\*nS(w' - ^^^nlX^) (X W \^S(W - r ^)\ 75*n|0>- 



(15) 



The computation of J 1 ^ is straightforward. At lowest order in a s , the momentum fractions 
w and w' should be the same because there is no collinear gluon emission to change the final 
momentum fraction. In this case J 1 ^ is given by 



with^ ( %,;4) 



J* (w,v/,px) = 5(w-w')K*' m (w,p 2 x ), 
n(l — w). 



(16) 



Putting Eq. (15) into Eq. (14), we finally obtain the factorization formula as 

F^(Q 2 ,x) = ^(^Y C du'dw'H*(u',w',Q 2 )Mu) f dudwH(u,w,Q 2 )<P^u) (17) 
2 , K y Q / Jo Jo 

xJ«(w,w',Q 2 (l-x)), 



with p\ = Q 2 (l — x). As seen from Eq. ( 15 ), J^ is the quantity of order 1, but it can include 
the logarithm of ln(Q 2 (l — x)/fi 2 ) at higher orders in a s . Therefore F£ is power counted 
as f 2 /Q 2 ~ A 2 /Q 2 ~ (1 — x) 2 because all the other quantities are of order 1. This power 
counting is consistent with the result in Table [I] 



B. Contribution from hard-collinear gluon exchange 

For hard-collinear gluon exchange, we employ two-step matching procedure QCD— > 
SCETi — y SCETn by integrating out hard {p\ ~ Q 2 ) and hard-collinear (p1 c ~ QA) degrees 
of freedom in turn. In SCETi we do not distinguish the hard-collinear and the collinear 
fields allowing the fluctuations of QA, while we keep only the collinear and soft fields with 
the fluctuations of A 2 in SCETn after integrating out the hard-collinear fields. 

At tree level the electromagnetic current operator J M = qj^q in the full theory can be 



expanded in powers of A = JA/Q in SCET: as 



2i )V _ . 2?; M — / 

- —y n Y^Y n $ nl _m n - *^ ± y_ty A + (\ 2 ), (18) 

where B% = [WliD£W n ], B£ = [WliD^Wn] are the gauge-invariant collinear gauge fields, 



and the derivative operators act only inside the bracket. The first term in Eq. (18) is the 
leading current operator, the remaining operators are of order A. 

Now we consider the hard-collinear gluon exchange between the electromagnetic current 
and the spectator quark. The spectator interaction is described by the following soft-collinear 
Lagrangian [TJ)J [16] 

£« = ^nftYiq. + h.c, (19) 

4 2 c ») = * n |n • B n Y^q 9 + h.c, (20) 

4? = *n| ^^M?. + h.C, (21) 

where the superscripts (i) in C sc denote the \ l suppression compared to the leading SCET 
Lagrangian. 

The contribution of the hard-collinear gluon exchange is described in terms of the time- 
ordered products of the electromagnetic current and the soft-collinear Lagrangians in SCETj. 
However when we go down to SCETn after integrating out the hard-collinear degrees of 
freedom, the power counting changes accordingly. The collinear momentum in SCETn scales 
as p M = (n ■ p,p±,n ■ p) = Q(l,r],7] 2 ) with 77 = A 2 . The power counting of the n-collinear 
fields and their derivatives, V± and n ■ V changes (^ n ,'P±,'n • V) ~ (A, A, A 2 ) — > (t),T],T] 2 ) 
when matched onto SCETn. 

This fact implies that the final power-counting in SCETn can be different from the power 
counting in SCETi. An example is a hard-collinear gluon exchange in SCETi from the 



time-ordered product of the leading electromagnetic current in Eq. (18) and Cy}, with the 
leading collinear Lagrangian Cf\ C^ contains an operator with two -D^'s, from which one 
hard-collinear gluon is contracted with C^J, and V± is selected from another. The resultant 
operator in SCETn has a V±, acting on the external n-collinear fielcjjand it is suppressed 



2 The derivative operator, V± docs not vanish unless it returns a total transverse momentum of the pion. 
The nonvanishing V\_ contributes to twist-3 LCDAs if we take the matrix elements for the pion [TS] . 



10 



by A. Therefore this contribution is eventually power-counted as the same order as the 
operators from subleading time-ordered products, which include neither V± nor n ■ V. As 
a result, we include all the subleading contributions in the time-ordered products in SCETi 
and the spectator contributions are given as 

n = ijd*x T{jj°^(0), £«(*)}, n = ijd'x T{J^(0),jrW (*)}, (22) 

n = -Jd'xd'y T{J^(0),CM(x),£M(y)}, I? = ijd'x T{j£*(0), /£>(*)}, 

where JJ 0)At = ^ n Y^Y n ^ n and J { L 1]>1 = -(2v fl /Q)^ n Y^Y n $ nl _m w which are the first and 
the fourth operators in Eq. (18). CJ£ is the subleading collinear Lagrangian given by [T7| 



£« = * n YtipJY n ^WlipiW n * n + h.c. (23) 



The fifth operator in Eq. ( 18 ) describes the interaction of n-collinear particles, that is, the 
jet function, and it has been considered to give a dominant contribution to the longitudinal 
structure function F L , while its overall contribution is suppressed compared to Fi [7HH]. But 
it turns out that J L is another source for Fl, and this contribution to Fl is comparable to 
the suppressed F\ near the endpoint region in the power counting from the above analysis. 



Our approach to the leading contribution in Eq. (22) is similar to the analysis for the 
heavy-to-light current for B — > it or K transition in semileptonic B decays [HI [20] , where 
the leading and subleading current operators involving a collinear gluon give comparable 
contributions in the power counting of l/m&. The leading current obeys the heavy-to-light 
spin symmetry [21], but the matrix element for the time-ordered products is nonf actor izable. 
It also has an endpoint divergence [19] or large ambiguities [20] . The remedy for this problem 
is to absorb the nonfactorizable contributions to the form factor. However the contribution 
of subleading currents violates the spin symmetry, but it is factorizable. In DIS with an 
initial pion, j\_ and J L have also different spin structures. In a similar manner the time- 
ordered products of j\_ have endpoint divergences if we take LCDAs for the pion, and they 
are absorbed into the nonperturbative hadronic matrix element, while the time-ordered 
products of J L give factorizable contributions, and are free of endpoint divergence. 

Including the radiative corrections, the relevant electromagnetic current operators can be 
written as 

4 1)M = ~ Jdu C L (Q,u,fi)y n Y n Y n ft n± 5(u- ^)*n, (24) 

where Ci(Q, fi) and Cl(Q, u, fi) are the Wilson coefficients. The Wilson coefficient Ci(Q, //) 
has been computed to one loop [6] . Note that Cl{Q, u, //) depends on the momentum fraction 
u of the incoming quark because J L is a three-particle operator. The anomalous dimension 
of Cl is given by Eq. (C8) in Ref. [7] to one loop. The renormalization group behavior of 
J}_ and J L is different since they are not a reparameterization-invariant combination [22J. 
Schematically both the structure functions F\ and Fl can be written as 

Fi~Hi-J n ®K u (25) 

Ff c ~H L ®Jn®J L ®S L ® $ L , (26) 

11 



where Hi = |Ci(<5)| 2 , Hl = Cl(Q, u)Cl(Q, v) are the hard factors, ® denotes an appropriate 
convolution. K\ is the hadronic matrix element of collinear and soft operators, which come 
from Tf, T% and T£. The contributions from Tf , T^ and T% can be written in the form 
Ji ® Si (g) $1, but they contain endpoint divergences. On the other hand, Fjj c from T£ is 
factorized. We put the superscript l HC on F L to distinguish it from F^, the hard gluon 
contribution to Ft- There are two kinds of jet functions Jn(x,fi), and Ji,l(x,^), obtained 
by integrating out the degrees of freedom of order p 2 ~ QA in the n and n directions 
respectively. Physically Jn(x,fi) describes the final states, while i7i,l(w, /-0 describes the 
initial states of collinear particles. Si t L are the soft functions including soft Wilson lines and 
soft spectator quarks, and $1^ are the LCDAs squared of the pion. 

In the conventional approach without including the spectator quark, the structure func- 
tion F\ can be cast into the following factorized form 

> n-p H - h 

n-p H 

where f q / n is the standard PDF obtained from the matrix element of a gauge-invariant 
collinear quark bilinear operator, 



F 1 (Q 2 ,x) = Hi(Q 2 , /i) J dlJ n (Q(l - x) - I, jf )/„/„( - 'J H l ), (27) 



f q/H (y) = (H\^ n ^6(yn -p H -n- V)* n \H), (28) 



and Hi and Jn are given in Eq. (25). The PDF can be additionally factorized into the soft 
and n-collinear parts, the combination of which recovers the renormalization behavior of the 
PDF [3 [23]. 



The factorization formula, Eq. (27), holds even when the spectator contributions are 
included. It can be achieved if we generalize the definition of f q / n with the spectator contri- 
bution K\. That is justified because the spin structure is the same for both contributions, 
and furthermore the renormalization behavior is also the same. Note that the structure func- 



tion F\ is scale independent and the remaining parts in both the expressions of Eqs. (25) 



and (27) are the same, therefore the renormalization group behaviors of f q u and Ki are also 
the same. In other words the spectator quark contributions involved in Ki are described by 
C sc in SCETi, which is scale independent and does not affect the renormalization behavior. 
Therefore we can safely generalize f q / n to Ki without inducing additional complications, 
and Ki or the PDF can be treated as a nonperturbative function to be determined from 
experimental data. As a result the definition of the standard PDF is still applicable near 
the endpoint region, but it holds up to SCETi. If we go further and employ the two-step 
matching, when the PDF is matched onto SCETn including the spectator contributions, it 
has more complicated substructure involving the lightcone distribution amplitudes of the 
initial hadrons. Note that Ki or the PDF can be dependent on the scattering processes es- 
pecially due to the difference of the soft functions in each scattering process. Theoretically 
the two-step matching result is more explicit, but it is not economical to express a nonper- 
turbative quantity Ki in terms of the convolutions of other nonperturbative quantities such 
as the LCDAs. 

For F^ c , we introduce a new nonperturbative function fi to cover Jl <g> Sl ® $l in 



Eq. (26). Note that the renormalization behavior of fi is different from f q / n because Hi 
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FIG. 2. Examples of the time-ordered products for the hard-collinear gluon exchange. Diagram 



(a) denotes T^ and (b) describes T^ in Eq. (22). The solid lines are collinear fermions, the dotted 
line denotes an usoft quark, and the wiggly line with a solid line is an n-hard-collinear gluon with 
P 2 ~ QA. 

and H L have different anomalous dimensions. So ft, is not related to f q /^, and it is a new 
contribution to the PDF near the endpoint region. As we notice in the case of heavy-to-light 
transition in B decays, fi can be factorized without the endpoint divergence. 

Now we consider the factorization proof for Fl in detail in order to see how the spectator 
contributions can be treated in the inclusive scattering process. The first step is to compute 
the hard-collinear gluon exchange and construct a four-quark operator consisting of two 
incoming collinear quarks, an outgoing collinear quark and a soft quark. This four-quark 



operator with a soft quark can be obtained by the time-ordered product, T£ in Eq. (22). 
The corresponding Feynman diagram is shown in Fig. |2J (a). After integrating out the 
hard-collinear gluon at tree level, T£ is written as 



rpfJ, 

1 1 



8iTa, 



Q' 



du 



C L (u) r drj 



u 



-Jdv) 



(29) 



x * R Y£Y n T a 'ftV n • ^ nl i T a 5( V + n ■ id)Y^q s 



where u — 1 — u. And Jl is the jet function obtained by integrating out the hard-collinear 
gluon in the n direction, with the normalization Jl(v) = 1 + 0(a s ). At higher orders in a s , 
there can be a color singlet four-quark operator with the structure 1(8)1. Since the initial 
pion is a color singlet, we take the appropriate color projection. The matrix element of T£ 
is given by 



(X|7T|7T> 



IV 



^nCpots f w f du 



N Q 



-C L (U) 



[U 



a 



(30) 



x / ^J^iX^j^YJir, + in ■ d)Y^q m \0), 

T] 2 



where (f) w is the leading twist pion LCD A in Eq. (11) 



The contribution of the hard-collinear gluon exchange to Fl is obtained by replacing J^ 
by T% in Eq. (4), and the corresponding Feynman diagrams with J L before integrating 
out the hard-collinear gluon and with T£ are shown in Fig. M (a) and (b) respectively. The 
discontinuity of the Feynman diagrams in Fig. [3] yields the structure function. As a result 
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(a) 



(b) 



FIG. 3. The Feynman diagrams for the longitudinal structure function Fl with the hard-collinear 
gluon exchange, (a) The exchanged hard-collinear gluon in the n direction is shown, (b) The 
equivalent diagram to (a) in terms of the time-ordered product of T^ and T%. 



the factorized contribution to F^ is written as 

J n 



F^ c (Q 2 ,x,fi) = 2-^ I dudvT L (u,v, n)(f)^(u, fifaty, fi), 



Q 
where the kernel Ti(u,v,fi) is given by 

8-Ka^Cp Hl(Q, u, v 



(31) 



T L (u,v) 



N 2 uv J mu *W L - X) - L )Jj^f 

In the case where the active quark is an antiquark and the spectator is a quark, the con- 



J dlJk(Q(l -x)-l)J ^^Mv, t/)S l (1, rj, rf). (32) 



tribution is the same because of the charge symmetry. So we put the factor 2 in Eq. (31) 



reflecting this fact. Here the initial jet function is given by J^l(V:V') = Jl(v)Jl(v')- The 
soft function Sl{1,i],i]'), which consists of soft quarks and soft Wilson lines, is written as 

S L {1,7],7)') = (0\T[q s YJ(r ] , -n-i^)Y^Y n ]p(l + n-zd)T[YiY n 6(7 1 + n-td)Y^q s \\0), (33) 

where T denotes the anti-time ordering. The discontinuity of the soft quark propagator in 

2n5{k 2 )k; 



the soft function in Eq. (33) gives the factor 

d 4 k 



(34) 



(2tt) 
from which the soft function at leading order in a s is written as 

Sf\h V, rf) = I ^5(k 2 )2n ■ kS(rj' - n ■ k)6( V - n ■ k)5(l -n-k) = -^I5( V - rf). (35) 
J (Zir) 6 loir 6 

Note that the soft function is denned to be dimensionless. Definitely this is different from 
the soft function appearing in the conventional approach which consists of only soft Wilson 
lines. The presence of soft quarks gives a different soft function. And the final jet function 
Jn{n • Px„) with n • px w = Q(l — x) — n ■ p Xs {Pxs being the total momentum of the soft 
particles) is obtained from the relation 

d A px n f 

(2vr) 4 2' 



£<0|tf fl |X fl )<X fi |tf fl |0) 



A"n 



., „J(n-Px n 



(36) 
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In this notation, the jet function at tree-level is given by J(r}) = 2ir5(ri) and it has been 
computed to two-loop order 



Since f n is 0(A), F^ c in Eq. (31) is power-counted as rf ~ (1 — x) 2 as we expected in 
Table [11 From Eqs. (31) and (32), the new nonperturbative function fi reads 

8irf 2 a 2 a C 2 F </>„( u,/i)^ w (u ,//) rdrjdrf . , , 

h{u,v,fj,) = — ^ = = / r J L {ri,ri,fi)S L {l,ri,ri,fi). 37 

Because Hl and J dlJn are of order 1, /j, is also power-counted as order rf. The same 
reasoning leads to the fact f q / n ~ rf ~ (1 — a;) 2 because W^" ~ F\ ~ (1 — a;) 2 . 

The result can be extended to the case with an initial proton in a straightforward way, 
but it is definitely more complicated because there are more spectator quarks. If we consider 
the similar factorization formulae for the structure functions F^ L ~ FL 1L (x or <g>) Jn® fi,L 
we can do the power counting on the nonperturbative functions f\ t L- Because W^ u ~ F^l 
rf ~ (1 — a;) 5 as seen in Tableland Jn is identical with the one defined in Eq. (36), both 






the nonperturbative functions f\ and fi scale as (1 — x) 5 . From Ref. [23] we can read off 
the fitted scaling behavior of the PDF from DIS experiments. At the factorization scale 
fiF = 3 GeV the powers of (1 — x) in the PDFs read ~ 4 for the u valence quark and ~ 5 for 
the d valence quark. It is consistent with our results considering huge uncertainties coming 
from the radiative corrections and renormalization scaling evolution. 

When we consider the time-ordered products for the hard-collinear gluon exchange in the 
proton, the electromagnetic current should be expanded to order 0{\ 2 ) since all the spectator 
quarks interact with the active quark. For example we obtain the following operator at 0{\ 2 ) 
to give a leading contribution to the structure function 

J ( l ] " = ~ f MuiCifavrfBtfylftlSfa - |^)jfc] [*(«i " |^)*»]. ( 38 ) 

where C' ± (ui,U2,fi) is the Wilson coefficient given by l/(wi + M2) at tree level and P M is 
the momentum of the proton. Since this operator is proportional to 7^, the time-ordered 



product contributes to F\. The anomalous dimension is C[ is different from C\ in Eq. (24) 
and hence we need a new nonperturbative function different from the standard PDF f q / p , 
which is induced from the time-ordered products of the leading electromagnetic current 

t(0)m 
J_l_ . 

Even though Fl is comparable to F\ in the power counting of (1 — x), the precise estimate 
on the size should include the radiative corrections and the evolution of the operators. 
The dominant contribution to Fi comes from the part involving f q / p , which is regarded as 
totally nonperturbative because the factorized expression f q / p = J <g> S ® $ is not justified. 
When Q 2 is large, a s (Q 2 ) or a s (Q 2 (l — x)) are significantly small. In this case, factorizable 
parts can be considered to be higher order in a s compared to f q / p . If the factorizable 
contributions are dominant in Fl, the size of F\ can be larger than Fl, which needs to be 
verified from experiment. For nonleptonic B decays, a similar comparison can be performed 
using experimental data. [2Sj For an initial pion, we have seen that Fl is totally factorizable 
both for hard-collinear and for hard gluon exchanges. But for a proton, a more detailed 
analysis is necessary in order to compare the size of F-y and Fl in the endpoint region. 
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IV. DRELL-YAN PROCESS NEAR THE ENDPOINT 

Near the endpoint in DY process with 1 — X\ ~ 1 — x 2 ~ 77, the quantity r = Q 2 / s 

approaches 1 with the power counting 1 — r ~ 77, where Q 2 is the invariant mass of the 

lepton pair and s is the hadronic center-of-mass energy. The variables X\ and x 2 are defined 

as 

Q 2 Q 2 

Xl = 2iv? X2 = 2iv? (39) 

where P\ and P 2 are the momenta of incoming hadrons. In this limit, the final-state invariant 
mass becomes 

p 2 x = Q 2 (l + -----) -> Q 2 (l - Xl )(l - x 2 ) ~ A 2 , (40) 

v r X\ x 2 ' 

requiring that only soft particles be allowed in the final state. 

Since the phase space in this endpoint region is so small, it is not interesting experimen- 
tally, but it is a good region to study the factorization property theoretically. To increase 
the available phase space, we may think of relaxing the condition such that p\ ~ QA. This 
region can be reached if only one parton is near the endpoint region, say, 1 — x\ ~ 1 and 
1 — X2 ~ 77. However, since the scattering cross section is a convolution with respect to X\ 
and x 2 , it is also possible to have 1 — x±, 1 — x 2 ~ y/fj such that (1 — £i)(l — x 2 ) ~ r?, 
which corresponds to none of the endpoint region. Actually, the region both away from the 
endpoint region is favored compared to the case with one parton near the endpoint region 
due to the steep decrease of the PDF near the endpoint. This region might be interesting 
on its own, but we confine to the above endpoint region here. 

The differential scattering cross section is given by 

w ^ Q, wsi^, Far ' (41) 

where Fqy is the DY structure function, which is given by [II] 

F DY = - J ^-9(q°)6(q 2 - Q 2 ) J d± ze' 1 ** {H X H 2 \ f(z) J^)\H X H 2 ). (42) 



(2nf 

Here J M is an electromagnetic current and the momentum q is given by q = P\ + P 2 — px, 
where P\ y2 are the momenta of two incoming hadrons H\ )2 . In the power counting the 
product of the volume elements d A z and d 4 q yields order 1 irrespective of whether the region 
is near or away from the endpoint. Near the endpoint, the label momenta, when integrated 
over the momentum, yields a Kronecker delta, and the remaining d 4 q is of order A 4 , while the 
volume element is of order A -4 . And away from the endpoint, d A q ~ Q 4 , and d A z ~ 1/Q 4 - 
However, there is a delta function S(q 2 — Q 2 ), which is power counted asD~ l/(Q 2 r]) since 
the argument in the delta function is given by q 2 — Q 2 = s(l — t)(1 — 2p x s~ 1 / 2 /(l — r)) of 
order Q 2 rj in the center-of-mass frame. 

We first consider the power counting of the hard-collinear gluon exchange contributions, 
and some examples of the contributing Feynman diagrams are shown in Fig. |4} Since there 
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should be only soft particles in the final state, hard-collinear gluon exchange is needed for 
each final soft quark. Following the same power counting rule as in DIS, the hard-collinear 
contribution to the structure function, namely F§y is power counted as 



F%S -V-M-I-F 



DY 



1 



1 



QA V Q 4 A 2 * 

1 / 1 \2 

QA' VQ8 A 4 



(QA) 4 • (A 



3\2 5 



(1 



(43) 



for H X o = n, 



(Q 3 A 4 ) 2 • (A 3 ) 4 ~ V u ~ (1 - r) 11 for H h2 = p,p 



where F is the power-counting on the final soft quark states. 

In DY processes, there is no analog of final-state collinear particles in DIS. However, there 
is another interesting process to be considered as far as the power counting is concerned. 
That process is "double parton annihilation", in which two quark-antiquark pairs in the 
incoming hadrons are annihilated by exchanging the momentum of order Q 2 . This process 
is shown in Fig. |5j The spectator quark-antiquark pair with energy fractions of order 1 is 
annihilated and transfers the whole energy to one of the active quarks. This is similar to the 
case with double parton scattering, but there is the difference in the final states in double 
parton annihilation. Since the momentum transfer is of order Q 2 , the resultant operators 
become local. Furthermore, they are lower in powers of a s compared to the corresponding 
hard-collinear gluon exchanges. That is, these contributions are of order a 2 (Q 2 ) and a 4 (Q 2 ) 
at leading order for pions and (anti)protons respectively. Using the power counting analysis, 
the power counting of the structure function for initial pions and p, p is summarized as 



F^ Y ^V-M-I 

1 /In-' 



QA vo3, 

1 / 1 \2 



QA Vq 



(44) 



(QA) 



4 ^3 



1 - r) 3 for H lt2 = 71, 



(Q 3 A 



3 a ^2 _ v 7 _ ^ _ T y for H = p) - 





FIG. 4. Examples of the Feynman diagrams with the hard-collinear gluon exchanges for initial (a) 
pions and (b) proton and antiproton near the endpoint, in which the spectator quarks become soft. 
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The Feynman diagrams in Fig. [5] can be dressed with soft gluons for final-state soft particles, 
but careful analysis of power counting shows that emission of soft gluons does not alter the 
result of the power counting without soft gluons. One thing to note in Fig. [5] (a) is that the 
Feynman diagram, when rotated, is exactly the same as the one for the pion form factor. It 
is interesting that the pion form factor and the double parton annihilation in DY processes 
are related. 

The complication in DY processes near the endpoint lies in the fact that there exists no 
limiting process from the conventional approach, and the double parton annihilation is less 
suppressed both in powers of a s and 1 — r. Among the contributions from hard-collinear 
gluon exchange, there can be nonfactorizable contributions when we take the LCDA for the 
initial state. If these nonfactorizable contributions are dominant, we can arguably regard 
Fdy ~ (1 — t) 5 or (1 — r) 11 from hard-collinear exchange without additional suppression 
by multiple powers of a s (QA), as we considered on the estimate of the sizes of F\ and Fl 
in DIS. In that case, these contributions from hard-collinear gluon exchange can be numer- 
ically comparable to the hard gluon contributions resulting in double parton annihilation, 
treating a s (Q 2 ) ~ 1 — r. On the other hand, if the double parton annihilation is the major 
contribution near the endpoint region, its effect may be noticeable as we get away from the 
endpoint region. But note that the conventional leading contribution of order 1 becomes 
dominant away from the endpoint region, and all the contributions considered above become 
subleading and are negligible. In some region between the standard region and the endpoint 
region, the effect of the double parton annihilation may be noticeable. However, for precise 
estimate and comparison, a more detailed analysis is necessary. 

The conventional approaches neglecting the spectator contribution have proposed the 
following factorization formula [U El [271 [28] 

F DY = H DY (Q 2 ) f" -S DY (1 - z)f DY (-), (45) 

J T z ^ z ' 

where Hdy is the hard function of order 1, Sdy is the soft function consisting of the products 
of the soft Wilson lines, and fuy is the convolution of the parton distributions, which is 





(a) 

FIG. 5. Feynman diagrams for double parton annihilation in Drell-Yan processes with hard gluon 
exchange between initial-state (a) pions, and (b) proton and antiproton. The diagrams with the 
gluons attached to other fermions connected to ® are omitted. 
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given by 



M^ = /'fw.)W9- («) 



Since / cteiSW(l — z) in Eq. (45) is of order 1. The power counting of the structure function 
in the conventional approach can be performed through fey. Since f q /u scales as (1 — x) 2 
for the pion and (1 — x) 5 for the proton in DIS according to our analysis, F^y can be 



power-counted as (1 — r) 5 or (1 — r) 11 treating the range of the integration in Eq. (46) to be 



of order r\. Therefore the estimate of the size in the conventional approach seems to favor 



the power counting of the hard-collinear contribution in Eq. (43). However it is not clear 
whether we can justify the parameterization of the contributions from hard or hard-collinear 
gluon exchanges as the convolution of the PDFs. 



V. CONCLUSION 

High-energy scattering processes near the endpoint region are hard to analyze in experi- 
ment, but they offer an intriguing opportunity to disentangle the structure of factorization 
properties in QCD. In this paper, a power counting analysis is performed for the structure 
functions in DIS and in DY processes near the endpoint region to claim that there are new 
contributions from hard-collinear gluon exchanges to be included since they are comparable 
to the currently available leading contributions. 

An important feature in this analysis is to apply kinematic constraints of the endpoint 
region to classify the possible types of final-state particles, while the initial partons and 
hadrons are required to be on the mass shell p 2 ~ A 2 . The resonance region is defined as 
the final states with p 2 ~ A 2 , and the endpoint region is defined as those with p 2 ~ QA. 
According to this classification, DIS can have both the resonance region and the endpoint 
region, but DY processes have actually only the resonance region. 

The explicit factorization proof for hard-collinear gluon exchanges in DIS is interesting in 
itself, but it is also illuminating to compare this process with nonleptonic B decays into two 
light mesons. In the factorization proof for nonleptonic B decays [29], we have considered the 
contribution of the four-quark operators along with the spectator interactions since they are 
of the same order. In the spectator interaction, a hard-collinear gluon is exchanged between 
the four-quark operator and a spectator quark in a B meson, and the final-state particles 
become collinear to form mesons. The hard-collinear gluon exchange considered here in DIS 
is exactly the reverse process of this spectator interaction, in which the final-state collinear 
particles are the incoming partons, and the initial soft quark is the final soft particle, and 
the heavy b quark is replaced by the n-collinear final-state jet. The factorization property 
of various spectator interactions is similar in both cases, noting the difference between a 
heavy quark and an n-collinear particle. This, along with the comparison between the 
double parton annihilation in DY processes and the pion form factor, shows an interesting 
relationship among different processes. 

In DIS, the spectator interaction has the same power counting as the process with final 
n-collinear particles, hence it should be included to be consistent. However, in DY pro- 
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cesses, the spectator interaction exists, but it is suppressed compared to the double parton 
annihilation. This result is surprising, but here we have considered only the power counting 
of various contributions, and we have not tried to give numerical analysis of those since it 
belongs to a future work. The power counting analysis indicates the degree of suppression 
in powers of 1 — x or 1 — r, but the actual contributions also involve other parameters such 
as some powers of a s at different scales Q 2 , QA. Therefore a study on the precise estimate 
of various contributions is necessary to compare with experiment. 
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